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Abstract 

We study a possible dark matter candidate in the framework of a minimal anoma- 
lous U (1)' extension of the MSSM. It turns out that in a suitable decoupling limit the 
axino, which is present in the Stiickelberg multiplet, is the lightest supersymmetric 
particle (LSP). We compute the relic density of this particle including coannihila- 
tions with the next to lightest supersymmetric particle (NLSP) and with the next 
to next to lightest supersymmetric particle (NNLSP) which are assumed almost de- 
generate in mass. This assumption is needed in order to satisfy the stringent limits 
that the Wilkinson Microwave Anisotropy Probe (WMAP) puts on the relic density. 
We find that the WMAP constraints can be satisifed by different NLSP and NNLSP 
configurations as a function of the mass gap with the LSP. These results hold in the 
parameter space region where the model remains perturbative. 
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1 Introduction 



A great deal of work has been done recently to embed the standard model of particle 
physics (SM) into a brane construction [1, 2, 3, 4]. This research is part of the effort, 
initiated in [5], to build a fully realistic four dimensional vacuum out of string theory. 
While the original models were formulated in the framework of the heterotic string, the 
most recent efforts were formulated for type II strings in order to take advantage of the 
recent work on moduli stabilization using fluxes. Such brane constructions naturally lead 
to extra anomalous £7(1) 's in the four dimensional low energy theory and, in turn, to 
the presence of possible heavy Z 1 particles in the spectrum. These particles should be 
among the early findings of LHC and besides for the above cited models they are also 
a prediction of many other theoretical models of the unification of forces (see [6] for a 
recent review). It is then of some interest to know if these Z 1 particles contribute to the 
cancellation of the gauge anomaly in the way predicted from string theory or not. In 
[7] some of the present authors have studied a supersymmetric (SUSY) extension of the 
minimal supersymmetric standard model (MSSM) in which the anomaly is cancelled a 
la Green-Schwarz. The model is only string-inspired and is not the low-energy sector of 
some brane construction. The reason of this choice rests in our curiosity to explore the 
phenomenology of these models keeping a high degree of flexibility, while avoiding the 
intricacies and uncertainties connected with a string theory construction. For previous 
work along these lines we refer to [8, 9, 10, 11, 12, 13]. In this work we perform a 
consistency check of our model [7] by evaluating the thermal relic density of the model to 
compare it against the WMAP data. If the axino is the lightest supersymmetric particle 
(LSP), its relic density is too high with respect to the experimental data. This is why we 
favor a next to lightest supersymmetric particle (NLSP) with a mass close to the LSP. 
We show that an interesting scenario arise also for three particle coannihilation processes. 
In these two cases, the model is consistent with the experimental data. Moreover in the 
three particle case we find configurations in which the LSP and the NLSP do not need 
to be nearly degenerate in mass. In this case the mass gap between the two can be of 
the order of 20%. This is the plan of the paper: in Section 2 we describe our model. In 
Section 3 and 4 we find the LSP and study the axino interactions. Finally in Section 5 
we compute the relic density. Section 6 is a summary of our results. 



2 Model Setup 

In this section we briefly discuss our theoretical framework. We assume an extension of 
the MSSM with an additional abelian vector multiplet with arbitrary charges. The 
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anomalies are cancelled with the Green-Schwarz (GS) mechanism and with the General- 
ized Chern-Simons (GCS) terms. All the details can be found in [7]. All the MSSM fields 
are charged under the additional vector multiplet V^, with charges which are given in 
Table 1, where Qi, Li are the left handed quarks and leptons respectively while L^, D\, Ef 
are the right handed up and down quarks and the electrically charged leptons. The su- 
perscript c stands for charge conjugation. The index % — 1, 2, 3 denotes the three different 
families. H u d are the two Higgs scalars. 





SU(3) C 


SU(2) L 


U(l)y 


U(l)' 


Qi 


3 


2 


1/6 


Qq 




3 


1 


-2/3 


Que 


m 


3 


1 


1/3 


Qd< 


u 


1 


2 


-1/2 


Ql 


Ef 


1 


1 


1 


Qe c 


H u 


1 


2 


1/2 


QH U 


H d 


1 


2 


-1/2 


Q Hd 



Table 1: Charge assignment. 



The key feature of this model is the mechanism of anomaly cancellation. As it is well 
known, the MSSM is anomaly free. In our MSSM extension all the anomalies that involve 
only the 577(3), SU(2) and U(l)y factors vanish identically. However, triangles with 
U(l)' in the external legs in general are potentially anomalous. These anomalies are 5 

(1) 
(2) 
(3) 
(4) 
(5) 

All the remaining anomalies that involve U(l)'s vanish identically due to group theoretical 
arguments (see Chapter 22 of [14]). Consistency of the model is achieved by the contribu- 
tion of a Stiickelberg field S and its appropriate couplings to the anomalous U(l)'. The 
Stiickelberg lagrangian written in terms of superfields is [15] 
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5 We are working in an effective field theory framework and we ignore throughout the paper all the 
gravitational effects. In particular, we do not consider the gravitational anomalies which, however, could 
be canceled by the Green-Schwarz mechanism. 
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where the index a = 0, ... ,3 runs over the U(l)', U(l)y, SU{2) and SU(3) gauge groups 
respectively. The Stiickelberg multiplet is a chiral superfield 

S = s + iV20iP s + 2 F S - ida^dd^s + ^9 2 9a^d^ s ~ \o 2 9 2 Us (7) 
and transforms under U{1)' as 

V (o) y(o) + ,( A _ A t) 

S -> 5 - 4i 6 3 A (8) 

where 63 is a constant related to the Z' mass. We give the expansion of £5 in component 
fields only for the part that is needed in the following sections. Using the Wess-Zumino 
gauge we get 
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As it was pointed out in [8], the Stiickelberg mechanism is not enough to cancel all the 
anomalies. Mixed anomalies between anomalous and non-anomalous factors require an 
additional mechanism to ensure consistency of the model: non gauge invariant GCS terms 
must be added. In our case, the GCS terms have the form [16] 

Cgcs = -d4(vMD"VM-VMD°vW)wM + h.c.] e2 p + 
+d 5 [(VWD"VW - V<®D a vW) + h.c.] e2§2 + 



+d 6 Tr 



(10) 
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where n.a.c. refers to non abelian completion terms. The b constants in (6) and the d 
constants in (10) are fixed by the anomaly cancellation procedure (for details see [7]). 

For a symmetric distribution of the anomaly, we have 

A A{0) h w h A(1) tfh A(2) A(4) 
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It is worth noting that the GCS coefficients d 4i5i6 are fully determined in terms of the *4.'s 
by the gauge invariance, while the b^'s depend only on the free parameter b 3 , which is 
related to the mass of the anomalous U(l). 

The soft breaking sector of the model is given by 

= M - \ (M A<°)AW + fc.c.) - \ (^Ws + h.c)j (12) 
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where C^jf M is the usual soft susy breaking lagrangian while is the gaugino of the 
added U (1)' and ips is the axino. The axino soft mass term deserves some comment. Based 
on [17] one could naively expect that such a term is not allowed since the Stiickelberg 
multiplet is a chiral multiplet. This turns out not to be correct because the Stiickelberg 
multiplet couplings are not of the Yukawa type. In fact in our model both the axino and 
the axion couple only through GS interactions. It is worth noting that a mass term for 
the axion <fi is instead not allowed since it transforms non trivially under the anomalous 
U(l)' gauge transformation (8). 



3 Neutralino Sector 

Assuming the conservation of R-parity the LSP is a good weak interacting massive particle 
(WIMP) dark matter candidate. As in the MSSM the LSP is given by a linear combination 
of fields in the neutralino sector. The general form of the neutralino mass matrix is given 
in [7]. Written in the interaction eigenstate basis (ip°) T = (ips, A^, \f\ hP d , hP u ) 

it is a six-by-six matrix. From the point of view of the strength of the interactions the 
two extra states are not on the same footing with respect to the standard ones. The 
axino and the extra gaugino A^ dubbed primeino are in fact extremely weak interacting 
massive particle (XWIMP). Thus we are interested in situations in which the extremely 
weak sector is decoupled from the standard one and the LSP belongs to this sector. This 
can be achieved at tree level with the choice 
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The neutralino mass matrix M.^ becomes 
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where M s , M , Mi, M 2 are the soft masses coming from the soft breaking terms (12) 
while My(o) is given in eq. (17). It is worth noting that the D terms and kinetic mixing 
terms can be neglected in the tree-level computations of the eigenvalues and eigenstates. 
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Figure 1: (a) Gaugino-fermion-sfermion interaction vertex, (b) Axino-gaugino- vector 
interaction vertex. 

so that the axino is the LSP. This assumption is motivated by the interaction strengths 
of the two extra states. Since the axino interacts via the vertex shown in fig. lb (which 
is of order ~ g^gl/Mz') while the primeino interacts via the vertex in fig. la (which is of 
order ~ g a ), if we do not assume the two decoupling relations (13) and (15) the dominant 
contribution in the (co) annihilation processes is that of the primeino, which is of the type 
of a standard gaugino interaction. 

The decoupling in the neutralino sector implies also a decoupling in the gauge bosons 
sector. In our model there are two mechanisms that give mass to the gauge bosons: (i) 
the Stuckelberg mechanism and (ii) the Higgs mechanism. In this extension of the MSSM, 
the mass terms for the gauge fields for Qh u = are given by 



where M V (o) = 4& 3 g is the mass parameter for the anomalous U(l) and it is assumed to 
be in the TeV range. The lower dots denote the obvious terms under symmetrization. 
After diagonalization, we obtain the eigenstates 
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with M 2 being the gauge boson mass matrix 
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and the corresponding masses 
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Finally the rotation matrix from the hypercharge to the photon basis is 
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where i,j = 0, 1, 2. 



4 Axino Interactions 

The axino interactions can be read off from the interaction lagrangian (9). The relevant 
term, written in terms of four components Majorana spinors 6 , is given by: 

C = iyft&WXVrffkf, 7l(d M K (1) )*s + i^&WTsIVS rmv^s (25) 

where the b^f 1 coefficients are given in (11). The related interaction vertex Feynman rule 
is 

C^[r,Y\ik, (26) 
where is the momentum of the outgoing vector and the C^'s are 

C <2 » = ^Slbf (27) 

The factor (27) contains the parameters which are related to the anomalous U(l) 
(see eq. (11))- Therefore <C g a and the axino interactions will be extremely weak, 
being suppressed by an order of magnitude factor with respect to the weak interactions. 
At tree level there is only one type of annihilation diagram, represented in fig. 2. We 

6 The gamma matrices -f^ are in the Weyl representation. 
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Figure 2: Annihilation of two axinos into two gauge vectors via the exchange of a gaugino. 

denoted with p\ and p 2 the incoming momenta of the axinos while k\ and k 2 are the two 
outcoming momenta of the gauge bosons in the final state. We will concentrate on the 
case with two photons in the final state. In this case the result for the differential cross 
section is given by 



where u\ and uo 2 are the energies of the two outcoming photons. Each amplitude M.i is 
proportional to the related coefficient whose generic form is given in (27). The cross 
section (28), being extremely weak, cannot give a relic density in the WMAP preferred 
range. Thus we are forced to consider a scenario in which coannihilations between the 
axino and the NSLP became sizable. Several scenarios can be considered for the NLSP. 
We can split them into two major classes: one in which the NSLP is either a pure bino 
or a pure wino, and thus a coannihilation with a third MSSM particle is needed in order 
to recover the WMAP result, and one in which the NSLP is a generic MSSM neutralino 
with a non negligible bino and/or wino component. In both classes in order to have 
effective coannihilations the NSLP (and eventually the other MSSM particle involved in 
the coannihilation process) must be almost degenerate in mass. As a first example we 
consider a pure bino as the NLSP. The allowed coannihilation processes with the axino 
are those which involve an exchange of a photon or a Z$ in the intermediate state and 
with a SM fermion-antifermion pair, Higgses and W's in the final state. The diagram 
with the fermion-antifermion in the final state is sketched in fig. 3. The differential cross 
section in the center of mass frame has the following general form 



where s is the usual Mandelstam variable and pf ti is the spatial momentum of the outgoing 
(incoming) particles. On dimensional ground \A4\ 2 has at least a linear dependence on 
Pf and this implies that the dominant contribution comes from the diagram with the SM 
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Figure 3: Coannihilation of an axino and a bino into a // pair via the exchange of a 
photon or a Z . 

fermion-antifermion pair / and / in the final state: 

* s A (o) - // (30) 
The resulting differential cross section, computed in the center of mass frame, is 



% - £ c, ^S§rw) [K + Ml - + M ' Mz - + M - M%) (31) 

where the sum is extended to all the SM fermions (with mass m/) while c/ is a color 
factor. Details of the amplitude computation can be found in Appendix A. 



5 Axino Relic Density 

In this section we compute the relic density of the axino. The case of the axino as a 
cold dark matter candidate has been studied for the first time in [18]. As we said in the 
previous section we study two scenarios: the first in which the axino coannihilates with 
only one NLSP degenerate in mass (a generic MSSM neutralino), the second in which 
there is an additional supersymmetric particle (either a chargino or a stau) involved in 
the coannihilation process with the axino and the NLSP. 

Just to fix the notation we briefly review the relic density computation for N in- 
teracting species [19, 20, 21]. The Boltzmann equation for N particle species is given 
by: 

^ = -ZHn - ( a iJ v ij)( n i n j ~ n ? n 7) ( 32 ) 

where rii denotes the number density per unit of comoving volume of the species i = 
1, . . . , N (i = 1 refers to the LSP, i = 2 refers to the NSLP, and so on), n = J2i n ii H is 
the Hubble constant, a^- is the annihilation cross section between a species i and a species 
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j, Vij is the modulus of the relative velocity while n eq is the equilibrium number density 
of the species v. 

n? = g i (l + A i f 2 e- AiX ' (33) 

where gi are the internal degrees of freedom, A, = (m; — m\)jm\ while Xf = m\/T is the 
freeze-out temperature. Numerical computation gives Xf ~ 20 [20]. 

Eq. (32) can be rewritten in a useful way by defining the thermal average of the 
effective cross section 

N n eq n eq 

(W> = EMi>^^ ( 34 ) 

obtaining 

^ = -3iJn-(a e// t;)(n 2 -K'') 2 ) (35) 

where n e,J = £V n^ 9 . As a rule of thumb [22] a first order estimate of the relic density is 
given by 

n , 2 lO-^cm^- 1 

tt Y h 2 ~ ; ; 36 

We take into account the two cases N = 2 and N = 3: 

• N = 2 case. Assuming that the relative velocities are all equal = v we get: 
/ (2) v / x (viiv)/((T22v) + 2(a 12 v)/(a 22 v)Q + Q 2 



where Q = n e 2 /n eq . The first term in the numerator can be neglected because the 
axino annihilation cross section is suppressed by a factor (C ( - a - ) ) 4 with respect to the 
MSSM neutralino annihilations (see the previous section) and thus (auv) <C (o' 22 v). 
The second term involves the coannihilation cross section. Let us consider the 
case in which the NLSP is a generic MSSM neutralino (a linear combination of 
hrf, h^) with a non vanishing bino or wino components. As we saw 
in the previous section each amplitude is generically proportional to C^gi with 
% = 1,2. Without loss of generality we consider the diagram which involves the bino 
component VP^A^ — > // and a photon exchange in the intermediate channel, i.e 
the M 2 amplitude in (31). We get 

C 2 = (C (1) cosBwf = 2(b ( 2 1) ) 2 gtcos 2 e w (38) 

From the expression of the mixed U(l)' — U(l)y — U (l)y anomaly (see [7]) and from 
the eq. (11) we have the following relation 

,(i) _ 3(3Q Q + Q L ) 

h ~ 256A (39) 
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where b 3 = M z >/4:g . With the assumption M z > = 1 TeV as in [7] we finally get 



,2 



~ 5.76 x 10- u (3g Q Q + g Q L ) 2 GeV~ 2 (40) 



where e is the electric charge. We get similar expressions for the other three terms 
in (31). This result has to be compared to the typical weak cross section (a 2 2v) — 
10~ 9 GeV~ 2 . As long as the charges and the coupling constant of the extra U(l) 
satisfy the perturbative requirement 

So-(3Qq + Ql) 2 <16 (41) 

the following upper bound is satisfied: 
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in the case of a pure bino, while 



< l(T b (42) 
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< 10- 5 (43) 



in the case of a pure wino. The previous ratios can be arbitrarly small, thus we can 
not derive a lower bound, but they cannot vanish since this would correspond to a 
complete decoupling of the axino. In this case the lightest MSSM state is the LSP. 
Accordingly to eqs. (36), (37), (42) and (43) the relic density gets rescaled as [23] 
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We performed a random sampling of MSSM models in which the NLSP is a pure 
bino or a mixed bino-higgsino (the case of a pure wino falls back into the N = 3 
case due to the wino-chargino mass degeneracy) and we computed the relic density 
in presence of coannihilations using the DarkSUSY package [24]. These two situa- 
tions are easily realized in some corners of the mSUGRA parameter space. Thus 
in our scan we assumed this scenario in order to fix the pattern of the supersym- 
metry breaking parameters at weak scale. We emphasize here that this choice is 
completely arbitrary, and it is assumed only for simplicity, since in our model [7] 
the supersymmetry breaking mechanism is not specified. In the former case there 
is no model which satisfies the WMAP constraints [25]: 

0.0913 < tth 2 < 0.1285 (45) 

since the annihilation cross section of a pure bino is too low and the rescaling (44) is 
not enough to get the right relic density. In the latter case the higgsino component 
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Figure 4: Axino relic density in the case in which the NSLP is a linear combination 
bino-higgsino. Red (darker) points denote models which satisfy WMAP data. Left panel: 
A 2 = 1%. Right panel: A 2 = 5%. 



tends to increase the annihilation cross section and thus we find models which satisfy 
the WMAP constraints. The results are summarized in the fig. 4 for A 2 = 1% and 
A 2 = 5%. In order to fulfill the WMAP data (red (darker) points in the plot (4)) 
the axino mass must be in the range 50 GeV < Ms < 700 GeV in the limit A 2 — > 0, 
where the lowest bound is given by the current experimental constraints 



N = 3 case. This is the case in which there is a third MSSM particle almost 
degenerate in mass with the LSP and the NLSP. Typical situations of this kind 
arise when the NLSP and the next to next to lightest supersymmetric particle 
(NNLSP) are respectively the bino and the stau or the wino and the lighest chargino. 
Expanding in an explicit way all the terms in the sum (34) we get: 
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(a 31 v)nl q n e ^ + {(J 32 v)n 3 q n 2 q + (a 33 v)(n' : 
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[(a 22 v) (niy + 2{a 23 v)n eq n cq + (033^) (^ q ) 2 ] 



(46) 



(n ec i) 2 

where in the last line we have neglected the terms (cruv), (cr^v) and (ai 3 v) since 
these are the thermal averaged cross sections which involve the axino. By introduc- 
ing a new set of variables defined by 
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for i = 2, 3 
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where are the internal degrees of freedom of the particle species, Xf — m\/T and 
A, = (mi — mi) /mi, we obtain 

/ (3) v _ (^22^)^2 + 2((T 23 v)Q 2 Q 3 + (^33^)^3 ^gx 

cS - ( i + g 2 + g 3 )2 

Under the assumption (m 3 — m 2 )/mi <C Q3/Q2 — 92,1 92 we finally get 

("SM * 7 ; Ql ; -2 ^mssmv) (49) 



i + (i + S)Q 2 



where 

(o"mssm^) = (cr 22 v) + 2—(a 23 v) + ( — J (0-33^) (50) 

#2 \<72/ 

In order to compute the rescaling factor between the relic density of our model and 
the MSSM relic density we have to express omssm in terms of a two coannihilating 
species effective cross section. This is given by 
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since {m 3 — m 2 )/mi <C 1/xf and m 2 > mi then (m 3 — m 2 )/m 2 <C We remind 

the reader that the values of n 2 q , n eq and n eq are different with respect to those in 
the former case since now there are only two species in the thermal bath. We then 
find 
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and inserting back this relation into (49) we obtain 





(54) 



1 + 



i 1 + s) ^ 



The rescaling factor between the three and two particle species relic density is given 
by the following relation 



We performed a random sampling of MSSM models with bino-stau and wino- 
chargino coannihilations. The first situation is realized in some corners of the 
mSUGRA parameter space 7 while the second situation is naturally realized in 
anomaly mediated supersymmetry breaking scenarios. For each model we com- 
puted the relic density (Vth 2 )^ for the two coannihilating species with the Dark- 
SUSY package [24]. We finally computed (fi/i 2 ) (3) using (55). The bino-stau 
models which satisfy the WMAP constraints have an axino mass in the range 
100 GeV < M s < 350 GeV in the limit A 2 — > 0. As the mass gap increases the 
number of allowed models drastically decreases and eventually vanishes for A ~ 5%. 
In the wino-chargino case, models which satisfy the WMAP constraints are shown 
in fig. 5 for four reference values of A 2 . The space of parameters with A 2 < 5% 
and an axino mass Ms > 700 GeV is favored while as the mass gap increases lower 
axino masses become favored, e.g. 100 GeV < M s < 200 GeV (A 2 ~ 20%). 

6 Conclusions 

We studied a possible dark matter candidate in the framework of our minimal anomalous 
U(l)' extension of the MSSM [7]. In the decoupling limit (13) and under the assumption 
M 3> M s , M V (o) the axino turns out to be the LSP. Being an XWIMP the axino annihila- 
tion cross section is suppressed with respect to the typical weak interaction cross sections. 
This implies that in order to satisfy the WMAP constraints on the relic density we must 
have at least a NLSP almost degenerate in mass with the axino. We considered the case 
with two and three coannihilating particles and we found some configuration which sat- 
isfies the WMAP constraints. The results depend on the mass gap between the axino 
and the NLSP. In the exact degeneracy limit A 2 — > the allowed models have an axino 
mass in the range 50 GeV < M s < 700 GeV for the bino-higgsino coannihilation case 

7 or in the so called Constrained MSSM (CMSSM). 




(55) 
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Figure 5: Axino relic density in the case in which the NSLP is a wino while the NNLSP 
is the lightest chargino. Red (darker) points denote models which satisfy WMAP data. 
Upper left panel: A 2 = 20%. Upper right panel: A 2 = 10%. Lower left panel: A 2 = 5%. 
Lower right panel: A 2 = 1%. 



while 900 GeV < Ms < 2 TeV for the wino-chargino coannihilation case. When the mass 
gap is A 2 ~ 20% the allowed models are those with wino-chargino coannihilations and 
an axino mass of 100 GeV < Ms < 200 GeV. Finally let us comment on the differences 
between our scenario and that studied in the work [23]. In our framework the U(l)' does 
not arise from a hidden sector and thus all the MSSM fields can be charged under this 
extra abelian gauge group. This is the most relevant feature which could also be detected 
experimentally (see for example [6]). Moreover in our scenario the axino interactions are 
suppressed with respect to the weak interactions due to the GS couplings while in [23] 
the mechanism to suppress the couplings and give an XWIMP is provided by the kinetic 
mixing between the U(l)' and U(1)y- 
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A Amplitude for Ai + ips — * // 



In this Appendix we give some detail about the amplitude computation for the process 



up 



ft,- Z, ft/" - — iW Zn 

(56) 

where C 1 = C^cosOw, Cz = —C^sm9 w while k 2 = s is the momentum of the 
intermediate gauge boson. The corresponding square modulus is 



\M\ 2 = -64 



a \k 2 -M*J k 2 + ^-MlJ 



(57) 



with 



T v = mj(p XlPs ) + Mi M s 2m) + S(p f p f )m 2 f + (p fPf ) 



+ 



(pxiPs)(pfPf) + 



T 



-(PfPf) [(P\iPf)(PfPs) + (P\iPf)(PfPs) 
-2(P\iPf)(PfPs) - (PMPf)(PfPs) ~ ip\J>f){pfPs) - 2(p Xl Pf)(PfPs) 

(PxiPf)(PfPs) + (pXiP^iPfPs) m ) - M 1 M S [mj - (PfPf) 2 ] + 
-{PfPf) {P\iPf){PfPs) + (p^prfiPfPs) 



(58) 



where p Xl , ps, Pf and pj are the bino, axino and SM fermions 4-momenta respectively. 
Writing all the momenta in function of s and integrating over the solid angle we get 



s — Am 2 x 



a = c f (g 2 b?) 

- 2Mf + (AM 2 + s) M\ - QM s sM l - 2M| + s 2 + M 2 s 



(59) 



x 



x 



12tt (M| o - s) 2 s 5 / 2 ^Mf - 2 (MJ + s) M 2 + (M 2 - sf 

(2m 2 + s) (2cos9 w eq f (M| o - s) + sm6 w g Zo v f s^ + (sm6 w gz a f f s 2 (s - Am 2 ) 
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